Energy gap in graphene nanoribbons with structured external electric potentials 
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The electronic properties of graphene zig-zag nanoribbons with electrostatic potentials along the 
edges are investigated. Using the Dirac-fermion approach, we calculate the energy spectrum of an 
infinitely long nanoribbon of finite width w, terminated by Dirichlet boundary conditions in the 
transverse direction. We show that a structured external potential that acts within the edge regions 
of the ribbon, can induce a spectral gap and thus switches the nanoribbon from metallic to insulating 
behavior. The basic mechanism of this effect is the selective influence of the external potentials on 
the spinorial wavefunctions that are topological in nature and localized along the boundary of the 
graphene nanoribbon. Within this single particle description, the maximal obtainable energy gap 
is -Bmax oc irhv-p/w, i.e., « 0.12 eV for w —15 nm. The stability of the spectral gap against edge 
disorder and the effect of disorder on the two-terminal conductance is studied numerically within 
a tight-binding lattice model. We find that the energy gap persists as long as the applied external 
effective potential is larger than ~ 0.55 x W, where W is a measure of the disorder strength. We 
argue that there is a transport gap due to localization effects even in the absence of a spectral gap. 

PACS numbers: 73.22.Pr, 73.22.-f, 73.20.-r 



I. INTRODUCTION 

The continuing rise of graphene as a new and ex- 
ceptionally promising material that outperforms conven- 
tional metals and semiconductors has initiated an ongo- 
ing quest for new physical effects. This has also generated 
a multitude of exciting proposals for various technical ap- 
plications, which have recently been summarized in sev- 
eral reviews^ - — However, due to single layer graphene's 
gap-less energy structure, applications where consider- 
able on-off current ratios are indispensable arc limited at 
present. Proposals for the creation of a lattice anisotropy 
that would lift the sublattice symmetry,* 4 , or for the ap- 
plication of strain fields^ that also could open an energy 
gap have yet to be realized. At present, bilayer graphene 
or certain graphene arm-chair nanoribbons have to be 
utilized instead, if an energy gap is needed. In the latter 
case, narrow ribbons of special widths have to be fabri- 
cated so that, due to quantum confinement, an energy 
gap or at least a transport gap in disordered ribbons is 
formed. For ribbon widths below 30 nm, the spectral gap 
is larger than kT at room temperature^^ Recently, the 
effect of a transversal electric field on arm-chair ribbons 
has also been studied theoretically^ 

Within simple non-interacting particle descriptions, 
graphene zig-zag ribbons are metallic and the open- 
ing of a spectral gap is impeded by electronic edge 
states^— appearing in an energy range where for broad 
two-dimensional graphene sheets valence and conduc- 
tion bands touch. These edge states are sensitive 
to an Aharonov-Bohm flux and robust against edge 
reconstructions ! 13 ' 14 For interacting electrons, it was 
shown 1 ^ that zig-zag ribbons always have a gap due 
to edge magnetization and that a homogeneous exter- 
nal electric field applied across the ribbon causes a 
half-metallic stated By employing the ab initio pseu- 
dopotential density functional method^ the authors of 



Ref. Il6| studied the spin-resolved electronic structure 
of zig-zag graphene nanoribbons and the possibility of 
spin-polarized currents. The influence of electron trans- 
fer between the two edges on the half-metallicity of 
the nanoribbon subjected to the electric field was also 
investigated*^ Recently, a gapped magnetic ground state 
has been suggested to be due to an antiferromagnetic in- 
teredge superexchange. 19 

These advanced theories are extremely interesting for 
clean graphene zig-zag nanoribbons. However, there is 
still no conclusive direct experimental observation prov- 
ing the existence of an one-dimensional magnetic state. 
The latter may well be spoiled in reality by edge disor- 
der or adsorbent atoms. — Therefore, we try to clarify 
in this paper whether one can obtain a spectral gap al- 
ready within a single particle description. In order to 
substantiate the relevance of such a basic model for real 
graphene zig-zag nanoribbons, we investigate the influ- 
ence of edge disorder on both the spectral properties and 
the two-terminal electronic transport. In our work, we 
first investigate a very simple spinless continuum model 
that can be treated analytically and then we employ a 
tight-binding lattice model including edge disorder ef- 
fects, which we solve numerically. 

Based on the Dirac equation for two-dimensional elec- 
trons with Dirichlet boundary conditions imposed in the 
transverse direction, we study in section |TT] the influ- 
ence of an effective potential acting within narrow strips 
along the edge regions of the graphene zig-zag nanorib- 
bon. This set-up, as sketched in Fig. (TJ can be studied 
experimentally in a three-dimensional device by apply- 
ing voltages between the back-gate of a graphene zig-zag 
ribbon and two top gates, one at the right and the other 
one at the left edge, respectively. The results are given 
in section Hill For a perfectly antisymmetric electric po- 
tential (left side V and right side — V, see Fig. [1]), we 
find a symmetric spectral gap at the Dirac point, which 
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increases linearly with the applied voltage V. Increasing 
the potential further, the gap reaches a maximum value 
of ~ tt/w, where w is the width of the ribbon, and finally 
closes again for even larger V. The splitting of the edge 
state energies, however, still continues to rise with V. 

In order to check how these results are influenced 
by disorder, we study numerically a tight-binding lat- 
tice model in section IPVI calculate the two-terminal con- 
ductance, and investigate the influence of edge disorder 
on both the spectral and the transport gap. The for- 
mer survives for V being larger than 0.55 x W, where 
W 2 /3 is the 2nd moment of the distribution P(e) = 
1/(2W)8(W - |e|) that defines the disorder potentials 
assumed along the ribbon edges. A transport gap, how- 
ever, is still observable even for very strong disorder when 
the spectral gap is absent. 
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FIG. 1: Cross section of a graphene zig-zag nanoribbon of 
width w (dashed line) showing the effective potential V(y) 
created by external voltages applied between the back-gate 
and two top-gates. The top-gates extend in the x direction 
along the left and the right edge having widths i/l and w — yR, 
respectively. 



II. MODEL AND SOLUTION 

We study a zig-zag nanoribbon of graphene, infinitely 
extended in x-direction and with finite width w in y- 
direction (see Fig. [T]) using the Dirac-type equation ap- 
proximation. This Dirac-fermion approach usually de- 
scribes correctly the low-lying states around the neutral- 
ity point in graphene £i22r21 There are two inequivalent 
(Dirac) points at k = ±K (valleys) in the band struc- 
ture where valence band and conduction band touch. 
The wave function for wave vectors near K is writ- 
ten as (£(a, x), x)), where a and b denote the two 
sublattices of the graphene structure. Correspondingly, 
(ji(a,x),r)(b, x)) is the wave function in the other valley. 
Since the Hamiltonian does not mix the valleys near K 
and K, the Dirac equation separates and reads for the 
first valley (in the other valley, we have a corresponding 
equation with id x — > —id x ) 



E - V(y) 
HvF(id x + d y 



hv F (id x - d y ) 
E - V(y) 



Here, vp — 10 6 m/s is the Fermi velocity, E the en- 
ergy, and V the electrostatic potential depending only 
on the y-coordinate. In what follows, we set hvp = I 
but recover the units when showing our results in the 
figures. In the corresponding lattice model, the bound- 
aries at y = and y = w are considered to be of zig-zag 
type. Then, in a description in terms of the Dirac model, 
we have periodic boundary conditions in the ir-direction 
and Dirichlet boundary conditions in the y-direction such 
that £(a, x, y — 0) = and £(&, x,y = w) — 

In order to have a simple model that can be treated an- 
alytically, we consider a piecewise constant electrostatic 
potential V(y) = V for < y < y L , V(y) = for 
y L < y < y R , and V(y) = -V for y R < y < w (see 
Fig. []}. At the points yL and yn the potential jumps, 
giving rise to a singular electric field only at yL and yn 
within the two-dimensional graphene sheet pointing into 
the y direction and being zero otherwise. ±V denotes 



the strength of the potential of the respective potential 
steps. This special choice leads to a symmetric energy 
gap around the Dirac points at E = if i/l — w — yR. 

Due to the periodic boundary conditions in the IE- 
direction, it is convenient to take the Fourier-transforms 
£(a,x) = e tc i x £(a,q,y), £(6,x) = e lqx £(b,q,y) and make 
the following ansatz for the wave function £(a, q, y) in the 
three regions of y 



L + e ikLV + L_e- lkL y, 
£(o, q,y)={ M+e ik **y + M_e~^\ y L < y < y R 
R +e ik R (y-u>) + R_e- ik *(y- w ),y R < y 



where kL,R and ku are given by 



(2) 



k L ,R = V(ETV)2-q 2 , k M = ^E 2 -q 2 . (3) 

The roots are defined to be positive if E is large and the 
usual analytic continuation is taken otherwise. £,(b,q,y) 
is then obtained from the Dirac equation as 

C q* L L+e^y + q L L-e-^y, y < y L 

Z(b, q,y) = { q* M M + e ik **y + q M M^ k ^ ,y L < y < y R 
{ q^R+e 1 ^-^ + q R R-e- lkR{ y- w \y R < y 

(4) 

with 

Q + ikL,R q + ikm , r , 

qL,R = ^_ TA , qM = — ^ — • (5) 
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E 



Here, q* M is given by q M \k M ^-k M even for imaginary k M 
and we have q* M qu = 1< The same applies to qL,R- 

The six amplitudes L±, M±, and R± follow from the 
normalization, from the boundary conditions 

C(o,g,0) = L+ + L_ = 

ZQ>,q,w) = q R R+ + q R R- = 0, (6) 

and from the matching conditions for £(a, g, y) and 
£(b, q, y) at y L ,R 

e L L + + e* L L_ = M'+M'_ 
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e M M' 



q* L e L L + + q L e* L L^ 
q* M e M M' + + q M e* M M'_ 



M'_ = e R R+ + e R R- 
= q* M M' + +q M M'_ 

= q*R e R R + + qr^rR- (7) 



with the abbreviations M' ± 



Jk L y L 



M±e ±ik M VL and 
ik M (VR-yL) „„ _ „ik R (w-y R ) 



e L = e , cm = e"" mKan aiJ ', e R = e 

(8) 

The "complex conjugates" e* L , e* M , and e* R are defined 
as for the qL, qM, and q R . Normalization of the wave 
function demands a non-trivial solution of Eqs. ([51 and 
that determines the energy. Furthermore, there are non- 
trivial solutions fci = leading to cl = qL = = 
— L_,M± = 0,R± = and similarly for k R = 0, k M — 
0. These result in a zero wave function and have to be 
excluded. 

After a straightforward calculation of the determinant 
of the 6x6 system, we get the condition for the eigenen- 
ergies in the form f(E; q, V) — with a real /: 



f(E;q,V) = 
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This function is even in V and for a symmetric arrange- 
ment, y L = w - t/fl, we have f(-E;q,V) = f(E;q,V). 
In the other valley, the eigenvalues are determined by 
f(E; —q, V) — 0. We denote the energy eigenvalues 
resulting from f(E;q,V) by E s<n (q), s = ±1, n = 
0, 1, . . . , oo. 

In the absence of an external electric potential (V — 
0), those wavefunctions corresponding to the eigenvalues 
close to E = 0, i.e., E Sj0 (q — > oo) with s = ±1, are 
localized along the edges^ With M± = L± and n = 
\/ q 2 — E 2 , one gets 

£(a,q,0 <y <w) = -2M+swh(ity) (10) 
£(b,q,Q < y < w) = s2M + sinh(n(w - y)). (11) 

Therefore, when the edge atoms at the left hand side 
are part of the b-sublattice and the edge atoms at the 
right hand side part of the a-sublattice, the wavefunction 
components of the b-sublattice are concentrated on the 
left (y = 0) for both energies E Sv0 (q), while we find the 
opposite for the components of the a-sublattice, which 
are concentrated at the right edge at y = w. The edge 
states' width depends on the imaginary momentum k. 
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FIG. 2: (Color online) The energy spectrum E a ,n(q;V) of 
a zig-zag graphene nanoribbon of width w — 15 nm in the 
absence of an external potential, V = 0. 



III. RESULTS AND DISCUSSION 

Now we turn to the discussion of the energy spec- 
trum. In Figs. [2H5J we show a numerical evaluation of 
f(E; q, V) — 0. In graphene zig-zag nanoribbons, the 
bulk gap closes due to surface state o 10 ' 11 at the edges at 
y = and y = w. This is shown in Fig [21 where part of 
the energy spectrum \E s , n (q)\ < 1.0 eV as obtained from 
© is plotted for one valley and V = 0. If we then apply 
a potential ±V that affects the eigenstates located at the 
the zig-zag edges as described above, a gap opens pro- 
portional to V until it reaches a maximum value deter- 
mined by the width w in y-direction. Here, we study the 
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FIG. 3: (Color online) The energy spectrum E Sl7l (q;V) of 
a zig-zag graphene nanoribbon of width w — 15 nm in the 
presence of external electric potentials V = ±0.0403 eV. The 
width of the left and right strips (jjl — w — yn) where the 
potentials are applied is 2w/5. 
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FIG. 4: (Color online) The energy spectrum E s ,n{q;V) of a 
zig-zag graphene nanoribbon for external electric potentials 
V — ±0.0978 eV. The widths of ribbon and potential strips 
are as in Fig. [3] 



symmetric case where — w — dr. For a finite electric 
potential ±V, the magnitude of the spectral gap depends 
on the width of the electrodes, ul — w — y^, as long as yL 
is smaller than the effective width w e wf 1 of the edge 
state. Therefore, a much stronger V must be applied to 
achieve the same spectral gap if j/z, < w e . The energy gap 
increases with yL and saturates at yL ^ w e for which the 
energy spectrum becomes independent of yL- The latter 
situation is realized in Figs. 02 IH andJSJ where the steps 
that define the width of the potential strips are chosen 
to be at yL — 0.4 w and yn = 0.6 w, respectively. If we 
continue to increase V, the gap closes again (see Fig. [5]) 
although the edge states, which can always be identified 
by their nearly dispersionless eigenenergies, move further 
apart because they are strongly affected by the external 
potential. With increasing V, the transition point where 
the bulk state transform into edge states moves to larger 
<?• 

The energy spectrum at q — can be found for arbi- 
trary potential widths yL and w — yn by putting in (|9]) 
Ql = Qm =qR = i- We get 

E s , n {q = 0) = s— (n+ -) - V . (12) 

w 2 w 

For n > 1 , the minima of the electron subbands and the 
maxima of the hole subbands appear not at q = but at 
q ~ 1/w. Also, for 2V < irhvp/w and increasing q, the 
electronic states at E S fi(q) start to become almost disper- 
sionless surface states close to the momentum q ~ 1/w. 
The corresponding eigenenergies can be calculated in the 
range 1/w < q < 1/ao, where ao is the lattice constant 
of the underlying lattice model. A careful evaluation of 
© in the limit q —> oo yields for the upper (s = 1) and 
lower (s = —1) state 
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FIG. 5: (Color online) The energy spectrum E 3 ,n(q;V) of 
a zig-zag graphene nanoribbon with electric potentials V = 
±0.1725 eV. A closer inspection of the spectrum shows that 
the curves do not cross. The widths of ribbon and potential 
strips are as in Fig. [3] 



( s\V\(l - e~ 2qVL ) ■ sV > 
E., (l/w < q < 1/ao) oc | e . 2q(w - yR)) J gV < Q 

(13) 

The appertaining density of states can also be esti- 
mated in this limit. It drops from large values at ±V 
down to zero at the Dirac point. For V = 0, however, we 
get E s ^(q) = s2qe~ qw and so the corresponding density 
of states behaves in the limit of \Ew\ -4 as ~ \Ew\~ 1 . 

Please note that in the Dirac approximation, the al- 
most dispersionless states at E±i (q) seem to be true 
solutions for all q — > oo. However, if one notices that 
this model originates from a lattice model, one sees that 
the nearly dispersionless states merely connect the two 
valleys at k = ±K and k = — K. That means, the Dirac 
model needs to be supplied with a cut-off (<7 max ~ 1/ao) 
in order to properly describe the asymptotic (E ~ 0) 
region of the tight-binding (TB) model. 



IV. TRANSPORT, INFLUENCE OF DISORDER 
AND CONCLUSIONS 

To support our finding of an induced spectral gap 
obtained within the continuum Dirac-model and to see 
its influence on the transport properties, we calculate 
numerically the two-terminal conductance g(E) of nar- 
row graphene zig-zag ribbons applying a transfer-matrix 
method within a TB-lattice model^ Here, semi-infinite 
leads are attached to both ends of the finite nanoribbon 
and the electric conductance is defined as usual via the 
transmission through the entire system (see Ref. [26| and 
references therein for details). As an example, the log- 
arithm of the energy dependent g(E) is shown in Fig. [S] 
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FIG. 6: (Color online) The logarithm of the two-terminal con- 
ductance g(E) vs. energy E of a perfect (blue curve) narrow 
zig-zag graphene lattice of 15 nm width in the presence an 
effective electric potential V/t — ±0.006 applied along the 
edges showing a transport gap AE — 2V. The red curve re- 
flects the influence of additional edge disorder W/t = 0.035 
(see text below). 



FIG. 7: (Color online) The transport gap AE vs. external 
potential V obtained numerically from the energy dependent 
two-terminal conductance g(E) of narrow zig-zag graphene 
lattices having widths 10 nm (▼), 15 nm (■), and 22 nm (•), 
respectively. The straight line is AE = 2V and t m 2.7 eV is 
the tight-binding hopping energy. 



for a clean 15 nm narrow nanoribbon with an electric 
potential V/t — ±0.006 applied along the edges lead- 
ing to a transport gap around the Dirac point. Here, 
t w 2.7 eV is the nearest-neighbor hopping term in 
graphene. The conductance exhibits sharp resonances 
with maxima g(E) ~ e 2 /h and decreases down to very 
small values ~ 10~ 10 e 2 /h for energies between ±0.006 t. 
The latter is due to tunneling and depends on the length 
of the sample. 

Fig. [7] shows the dependence of the induced transport 
gap AE on the effective external potential V, and that 
nicely confirms the linear behavior of our analytical re- 
sults for small V. The transport gap data plotted vs. 
potential in Fig ([7]) can be re-scaled by the respective 
ribbon widths. The outcome of this procedure is, within 
the given uncertainty of the data, a single fitting curve for 
all data (not shown). This result can be understood from 
an evaluation of ((9]) where one finds for yi, = w — i/r that 
the maximal gap appears always at the Dirac point q = 0. 
Therefore, Eq. (fT2j) can be used to see that 2Ei } qw = it 
which together with the observed 1/2 AE — V relation 
gives the scaling behavior mentioned above. The maxi- 
mal transport gap observed in the lattice model agrees 
within the numerical uncertainty with the spectral gap 
of both the finite lattice model and the continuum model 
where L — > oc was assumed. Recovering the units, we get 
from ([12]) and hv F = 3/2tc, where c = 1.42 x 10~ 10 m is 
the carbon-carbon bond length, that the maximal spec- 
tral gap (for ul — w — yn) follows the relation 

7T 3 7T 

E mM w = —nvF = -— tc, (14) 
w 2w 

leading to E max w 0.12 eV for W = 15 nm. We also 
find in our numerical work that replacing the assumed 



piece-wise constant effective electric potentials by more 
realistic smooth potential steps does not modify the re- 
sults. 

Finally, to check the robustness of this proposed gap 
opening mechanism against edge disorder, which may 
arise, e.g., through edge passivation by randomly placed 
hydrogen atoms that is known to stabilize the edges of 
pristine zig-zag nanoribbons considerably— we apply a 
random disorder potential along the border of the rib- 
bon. The eigenvalues are obtained by standard diago- 
nalization of the Hamilton matrix for L x x L y graphene 
zig-zag ribbons described by the TB-Hamiltonian defined 
on a bricklayer lattice 2 ^ with sites r and nearest neighbor- 
distance a 

H = 22e F ctc r -t ^2 c t c r>> ( 15 ) 

<r#r'> 

where (r ^ r') are pairs of those neighboring sites that 
are mutual connected on the bricklayer. The disorder 
potentials e r are uncorrelated random numbers that are 
non-zero only at the outer sites (different sublattice on 
left and right edge) along the zig-zag edges of the nanorib- 
bon and uniformly distributed between ±W, where W 
denotes the disorder strength. 

The resulting density of states (DOS) of a 15 nm zig- 
zag nanoribbon with L x /a — 720, L y /a — 72, and V/t — 
0.02, averaged over 1000 disorder realizations, shows a 
broadening of the DOS-peaks around E/t = ±0.02 that 
originate from the edge states (see Fig. [8]) . With increas- 
ing disorder potential strength W/t = 0.005, 0.01, 0.015, 
0.02, 0.025, 0.03, and 0.035, the spectral gap decreases 
linearly. This is seen in the inset of Fig. where the 
above results are shown together with additional data 
from a system of size L x /a — 960, L y /a — 48 (10 nm rib- 
bon), and V/t — 0.01. All data points collapse onto the 
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FIG. 8: (Color online) The density of states p(E) of graphene 
zig-zag nanoribbons (L x /a = 720, L v /a = 72) with ex- 
ternal edge potentials V/t = ±0.02 in the energy interval 
—0.025 < E/t < 0.025 around the Dirac point. One identifies 
the disorder broadening of the peaks at E/t — ±0.02, which 
belong to the eigenstates located along the ribbon edges, 
for disorder strengths W/t = 0.005 (black narrow peaks), 
0.01 (red), 0.015 (magenta), 0.02 (orange), 0.025 (blue), 0.03 
(green), and 0.035 (grey broadest peaks). The mean spectral 
gap Eg decreases linearly with increasing disorder strength W . 
This is shown in the inset (•) together with additional data 
for a system of size L x /a = 960, L y /a = 48 and V/t — 0.01 
(■) which all scale onto a single curve E g = 2V — 1.1 W almost 
until the gap closes. 



function E g — 2V — 1.1VF, which was also observed for 
other ribbon sizes having widths larger than 7 nm. There- 
fore, a gap should remain open in experiments when V 
is tuned to be larger than 0.55 x W, where W is usually 
fixed by the sample dependent intrinsic edge disorder. 
Here, the mean energy gap E g = l/N r Yl^i^t ~ ^i) 1S 
defined as the ensemble averaged difference between the 
smallest positive and the largest negative eigenvalue aver- 
aged over N r realizations. For general disorder with box- 
probability density-distributions P(e), we find a gap clos- 
ing relation E g = 2V — kT^ 2 , where T 2 = J M e 2 P(e)de 
is the second moment of the disorder distribution and 
k = 1.9 is an empirical constant. 

The influence of edge disorder on the logarithm of the 
two-terminal conductance, hi g(E, W) averaged over 100 
realizations, is shown in Fig.[B]for a finite lattice of width 
Ly/a — 72 and length L x /a — 720. Due to the ran- 



dom edge potentials, the sharp conductance resonances 
of the clean sample are smoothened out. Yet, a transport 
gap remains visible even in the case of strong disorder 
W = 0.035 1 when the spectral gap has completely van- 
ished but g(E) still drops six orders of magnitudes from 
about g ~ 0.1 e 2 /h at E = 0.008 1 to g = lO^V/Zi at 
E = 0. This means that the almost one-dimensional edge 
states of the clean sample become Anderson localized in 
the presence of sufficient edge disorder. This notion has 
been corroborated by an investigation of the respective 
eigenstates and by calculations of the length and disorder 
dependence of g(E). Previous studies have reached sim- 
ilar conclusions for nanoribbons with rough edgeSi 28 > 29 

In conclusion, we have shown that the application of 
external electric potentials, covering the area of the elec- 
tronic edge states that are located along the zig-zag edges 
of a graphene nanoribbon, can open a tunable spectral 
gap. Thus, one can convert the metallic behavior into a 
semiconducting one. For small potentials, the gap in- 
creases linearly with the potential strength, reaches a 
ribbon-width-dependent maximum Tthv-p/w (« 0.12 eV 
for w =15 nm) and closes again with further increasing 
electric potentials. The origin of this effect comes from 
the sensitivity of the spinorial edge states to electric po- 
tentials. Applying distinct external biases to the left and 
right edge state leads to a different shift of the almost dis- 
persionless edge energies as long as they are not pinned 
to the Fermi level. Using electric potentials of opposite 
sign causes the largest energy gap possible. The disorder 
effects, which may be due to atoms and molecules that 
saturate the dangling-bonds along the zig-zag edges in 
real samples, are found to reduce the spectral gap. The 
latter remains, however, finite as long as W < 1.82V, 
where W is a measure of the disorder strength and V the 
applied effective electric potential. For even larger dis- 
order strengths, a transport gap is still present allowing 
for reasonable on-off-rations for the electric current. Fu- 
ture experiments will show whether the present results 
of single particle physics are sufficient for the descrip- 
tion of a gap opening by external potentials in graphene 
zig-zag nanoribbons or if theories that emphasize edge 
magnetis m 15 ' 16 i 19 due to e-e interactions have to be ap- 
plied. 

Note added: During the review procedure, we became 
aware of a recent paper 30 by Bhowmick and Shenoy that 
addresses a spectral gap opening induced by external 5- 
like potentials placed along the edges of graphene zig- 
zag ribbons. This specific potential choice represents a 
special case contained in our model. 
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